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ABSTRACT 

Interstellar magnetic fields in elliptical galaxies have their origin in stellar fields that 
accompany normal mass loss from an evolving population of old stars. It is likely that these seed 
fields are amplified by interstellar turbulence which in turn is driven by stellar mass loss and 
supernova events. Since the local turbulent velocity is likely to exceed the global velocity of the 
interstellar cooling flow, magnetic fields are expected to be disordered and this is indicated by 
numerous observations. Further amplification occurs as these tangled fields are compressed in 
the inward moving cooling flow of the interstellar gas. Near the centers of galactic cooling flows, 
where the gas radiates away most of its thermal and gravitational energy, magnetic stresses are 
expected to dominate. We study here the time-dependent growth of interstellar magnetic fields 
in elliptical galaxies and demonstrate that fields similar in strength to those observed can be 
generated solely from natural galactic processes. 

Although interstellar turbulent velocities and therefore the equipartition field Be are larger 
near the galactic center, the interstellar field throughout elliptical galaxies can be determined 
by the outermost turbulent regions in the interstellar gas. This occurs whenever the increase 
in the field strength due to compression in the galactic cooling flow exceeds the rate of field 
amplification by local turbulence, i.e. when B > Bg. The magnitude of the interstellar field 
depends on the most distant radius in the galactic ISM at which the turbulent dynamo process 
can successfully have amplified stellar seed fields. Because of the long hydrodynamic flow times 
in galactic cooling flows, currently observed magnetic fields may result from periods of turbulent 
field amplification that occurred in the outer galaxy in the distant past. In our most optimistic 
turbulent dynamo models in which field reconnection is ignored, tangled interstellar magnetic 
fields of i? ~ I — 10 (r/IO kpc)“^-^^G are typical; these are consistent with observed fields. 
However, the sensitivity of the galactic magnetic field to poorly known historical turbulent 
conditions in the outer galaxy complicates estimates of field strengths that can be produced 
internally. Internal fields in ellipticals may also result from ancient galactic mergers or from 
shear turbulence introduced at the boundary between the interstellar gas and ambient cluster 
gas. 


Subject headings: galaxies: evolution - galaxies: magnetic field - galaxies: interstellar gas 
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1. INTRODUCTION 

It is likely that gas observed in the central regions of elliptical galaxies is dominated by magnetic 
stresses. The chaotic, non-equilibrium appearance of dusty clouds of cool gas apparent in HST observations 
in the cores of many elliptical galaxies (e.g. van Dokkum & Franx 1995) is often ascribed to the effects of 
a recent merger. However, it may be more likely that the disorganized appearance of these clouds reflects 
the well known inability of magnetic fields to reach quiescent configurations in differentially rotating, 
gravitationally bound fluids. In addition it is plausible that strong magnetic forces in the gas that cools near 
the centers of elliptical galaxies may influence the initial mass function or other aspects of star formation 
associated with galactic cooling flows. 

The interstellar medium in ellipticals and its magnetic field can be regarded as a result of cumulative 
stellar mass loss from an evolving population of old stars. Gas ejected from orbiting stars thermalizes a 
significant fraction of its kinetic energy and settles into a quasi-hydrostatic equilibrium with temperatures 
near the virial temperature of the galactic potential, Tyir ~ 10^ K. The hot gas then loses energy by emitting 
the thermal X-rays that are observed and slowly sinks in the galactic potential where it is compressed and 
maintained approximately isothermal. Although these subsonic inflows are commonly referred to as galactic 
“cooling flows,” the cooling occurs only in the central regions where radiative losses finally overwhelm 
compressional heating. 

The hot interstellar medium in elliptical galaxies operates like a magnetic distillery, first amplifying 
fields introduced in stellar ejecta then concentrating this magnetic energy in high density gas near the 
center of the galactic cooling flow where intense radiative cooling occurs. Small fields contained in winds 
and planetary nebulae expelled from red giants are thought to be amplified further by turbulent motions in 
the interstellar medium (ISM). Even in the presence of turbulent amplification magnetic stresses at large 
galactic radii are expected to be small compared to the pressure in the hot interstellar gas. However, the 
interstellar gas compresses as it slowly flows toward the galactic center and the field strength must also grow 
by flux conservation B ~ . The ratio of magnetic to (isothermal) gas pressure also increases with gas 

density, /p ~ p^/^. Finally, as gas cools near the central parts of the cooling flow, the thermal support 
of the gas drops precipitously and magnetic stresses are distilled out to support the remaining ionized gas 
against gravity and external gas pressure. 

Soker & Sarazin (1990) discussed the spherical, steady state evolution of completely disordered, 
frozen-in fields in cluster-scale cooling flows. They chose the field strength to be 1 pG at a radius of 50 
kpc. Within ~ 5 kpc in their models the magnetic pressure grew by compression in the cooling flow to be 
comparable to the gas pressure, but further amplification was not thought to be possible because of loss 
of field by reconnection. Since the cooling flow accelerates as the galactic center is approached, Soker & 
Sarazin noted that an initially random field would develop a large scale radial component. More recently 
Christodoulou & Sarazin (1996) computed the time-dependent evolution of large scale toroidal and axial 
fields in cluster cooling flows adapted to cylindrical geometry. Although B'^/Stt may be small relative to the 
local gas pressure, the origin of ~ IpG fields assumed by Soker & Sarazin and Christodoulou & Sarazin at 
the start of their calculations must be understood. Primordial field production prior to galaxy formation 
results in much smaller fields (^ 10“^® G) (Rees 1987; Lesch & Chiba 1995). But the high metallicity in gas 
in clusters of galaxies implies that much of this gas was enriched and expelled by Type H supernovae out of 
galaxies in which interstellar field growth mechanisms were already available. Thus cluster-scale magnetic 
fields probably have their origin in fields generated in the ISM of cluster galaxies. 

Our interest here is to explore the growth of interstellar magnetic fields in individual elliptical galaxies 
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in a self-consistent fashion beginning with small seed fields contained in stellar ejecta. We assume that 
stellar dynamos are efficient in producing fields in stellar envelopes. Stellar mass loss eventually transports 
these fields into the ISM where the held strength is greatly diminished by adiabatic expansion. However, 
recent theoretical work (Lesch & Bender 1990; Moss & Shukurov 1996) has suggested that some the kinetic 
and thermal energy that accompanies stellar mass loss processes in ellipticals can also drive turbulent 
motions in the ISM which in turn amplify the magnetic held to ~ 1 — 10 /iG. In this paper we shall assume 
with these authors that relatively small interstellar magnetic helds can grow by this turbulent dynamo. We 
describe how large scale helds may develop from random helds, but we suppose for simplicity that most of 
the held remains random on small spatial scales 0.1 kpc say). The growth of magnetic forces in the ISM 
may ultimately be limited by held reconnection which we allow for in an approximate way. 

Unfortunately, the physical processes that control the turbulent growth of galactic magnetic helds are 
poorly understood and controversial. In our Galaxy the mean held a — Q dynamo has been used to describe 
the amplihcation of small seed helds to equipartition with interstellar turbulence (Ruzmaikin, Shukurov, & 
Sokoloff 1988). Equipartition of magnetic and turbulent energy densities is in fact observed. Fields in spiral 
disks can grow both by turbulence and differential rotation. Enhancement of very small seed helds 10“® 
G) by the turbulent dynamo process is attractive because amplihcation by differential rotation in the 
Galactic disk is limited by the relatively few rotations (50 - 100 times) experienced by the Galaxy during its 
lifetime. However, several authors (Vainshtein & Gattaneo 1992; Kulsrud & Anderson 1992) have argued 
that the linear turbulent dynamo cannot amplify very small seed helds to equipartition values. Kulsrud 
& Anderson assert, for example, that the held strength can grow exponentially with time to equipartition 
only on small spatial scales where it is dissipated into heat while the mean held at intermediate and large 
spatial scales grows to saturation far below equipartition values. The (damping) reaction of the held on the 
background hydrodynamic turbulence is also considered. These authors claim that the failure of dynamo 
theory to account for helds of equipartition strength implies that Galactic helds must be primordial. But 
helds expected in proto-Galactic gas are likely to be much lower than those currently observed (Lesch 
& Ghiba 1995; Beck et al. 1996). Field (1995) has reviewed this apparent difficulty with the mean held 
dynamo theory and suggests that a more appropriate application of the Kulsrud-Anderson formalism would 
be to study the steady-state distribution of spectral held components in response to a constant input of 
turbulent energy on large scales. But an extension of the theory to include non-linear terms may also be 
necessary. In their recent comprehensive review of galactic helds Beck et al. (1996) note the apparent 
inadequacy of the linear theory in addressing turbulent growth; a non-linear turbulent dynamo theory 
would include inverse cascade of small scale held energy to larger scales as well as non-linear aspects of 
small scale dissipation. Finally, Beck et al. (1996) note that helds observed in the Galactic ISM are not 
concentrated to small spatial scales as the linear turbulent dynamo theory would predict. 

In view of these difficulties and possible shortcomings of mean held dynamo theory, we adopt here 
a simplihed model that allows turbulent growth of mean interstellar helds in elliptical galaxies. We are 
encouraged to take this deliberately optimistic approach by analogy with equipartition helds observed in 
the ISM of our own Galaxy and by the arguments of Moss & Shukurov (1996) and others that turbulent 
dynamos are likely to be effective in elliptical galaxies. By assuming the validity of turbulent dynamos, 
we can estimate the (maximum) held strengths expected in the ISM of elliptical galaxies from internal 
processes alone; no ad hoc initial galactic held is required. 

We have been motivated to study the evolution of magnetic helds in ellipticals by our recent work on 
the nature of cooling hows in slowly rotating ellipticals (Brighenti & Mathews 1996). We showed there that 
a small rotation of the stellar component, similar to those observed in most large ellipticals (~ 50 — 100 km 
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s~^), causes the interstellar gas to cool onto a very large disk comparable in size to the half-light radius. 
Most of the new gas reaching the disk arrives at its outer edge so that the cooling disk slowly grows with 
the age of the galaxy. As hot gas approaches the relatively high density regions near the outer edge of 
the cooling disk, its thermal X-ray image resembles a bright disk or torus in the equatorial plane. Near 
the equatorial plane the gas spins up to the circular velocity which is much larger (~ 400 — 500 km s“^) 
than that of the stellar system. We estimated that a small, dynamically insignificant magnetic held would 
amplify by a factor of > 5 in the differentially rotating flow as hot gas approaches the edge of the cooling 
disk. We speculated that the corresponding increase in magnetic forces, by > 25, might be sufficient to limit 
the maximum compression and X-ray visibility of the hot torus just beyond the edge of the cooling disk. 

If interstellar magnetic helds are sufficiently large, the X-ray images of rotating ellipticals may appear to 
be more circular than those predicted by Brighenti & Mathews (1966). This type of magnetic influence on 
the X-ray image may be necessary to account for the apparent circularity of Einstein and ROSAT images 
of several bright ellipticals in the Virgo cluster. Alternatively, the less battened X-ray appearance of these 
galaxies could simply result from a fortuitous orientation of galactic rotation axes along the line of sight. In 
any case, the growth and evolution of magnetic fields in spherical, non-rotating galaxies which we discuss 
here is a necessary hrst step toward a more complete study of held amplihcation and X-ray morphology in 
slowly rotating ellipticals. 

Our approach to held amplihcation is similar to that of Moss and Shukurov (1996) who describe the 
growth of magnetic helds in the turbulent interstellar environment energized by supernova explosions and 
the orbital energy of stellar ejecta. We extend the discussion of Moss and Shukurov to consider how the 
held is further increased (i) by compression as the cooling how approaches the galactic core and (ii) by the 
greater interstellar turbulence expected in past times when the stellar mass loss and supernova rates were 
both larger. Our calculation is necessarily approximate since many theoretical aspects of the turbulent 
dynamo process and held reconnection are based on dimensional arguments. Our models are exploratory, 
not dehnitive. 

We begin with a brief review of magnetic helds observed in elliptical galaxies then proceed to the 
description of our time-dependent models and the resulting magnetic held distributions. 


2. OBSERVATIONAL EVIDENCE FOR MAGNETIC FIELDS IN ELLIPTICALS 


The common signposts for the presence of magnetic helds in spiral galaxies - polarization of aligned 
dust grains, polarized emission from relativistic electrons and Zeeman splitting - are not available in 
ellipticals due in part to the relative absence of dust (Tsai & Mathews 1996a; b), the infrequency of 
supernova explosions and the absence of distributed cold HI gas. By far the best means of detecting helds 
in ellipticals is through the differential Faraday depolarization of radio sources by helds within the hot. 
X-ray emitting cooling how gas which lies along the line of sight. In most FR II-type double radio sources 
the jet that powers the radio lobes is visible at radio frequencies on only one side. Laing (1988) hrst noted 
that the radio polarization was systematically higher on the jet-side which is also the near side, assuming 
that radiation from the jet is made more visible by relativistic Doppler beaming. If near and far lobes are 
otherwise identical, the depolarization of the more distant lobe must arise in the galactic ISM, not within 
the lobes or in a shell of gas adjacent to the radio lobe. The body of evidence supporting this correlation 
has grown considerably, both for lobes in large ellipticals and in small clusters of galaxies (Carrington et 
al. 1988; Carrington, Conway, & Leahy 1991; Carrington & Conway 1991). The depolarization asymmetry 
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is particularly strong for projected lobe separations of galactic dimension, < 100 kpc. By comparing the 
relative Faraday depolarization of two lobes it is possible to determine the mean product of the electron 
density and net line of sight field strength to the center of the galaxy. A typical value is (ni?||) « 10“^ cm“^ 
/rG, implying i? ~ 1 — 10 /iG for plasma densities characteristic of the hot interstellar gas. These values are 
based on models in which the magnetic field is uniform in many small cells which are randomly oriented 
along the line of sight through the hot ISM (Burn 1966). Randomly directed fields are also implied by the 
patterns of differential radiofrequency Faraday rotation observed at high spatial resolution across the face 
of extended radio lobes or regions (e.g. Strom & Jagers 1988; Owen, Eilek, & Keel, 1990; Clarke, Burns, & 
Norman 1992). Magnetic fields even larger than those indicated above may be present if they are tangled 
on spatial scales much less than the radio resolution limit. Finally, Greenfield, Roberts, & Burke (1985) 
observed differential radio polarization in two gravitationally-lensed images of a distant quasar which they 
attribute to magnetic fields in the hot ISM of the intervening CD galaxy lens. 

We conclude that the hot gas in elliptical galaxies that participates in galactic scale cooling flows 
contains magnetic fields that are comparable to those in spirals and that a significant component of this 
field is spatially disordered. 


3. MAGNETIC COOLING FLOWS 


The equations that describe the evolution of hydromagnetic cooling flows in elliptical galaxies are the 
usual conservation equations modified with appropriate source terms: 


dp 

— + V • pv = ap*, 
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(3.3) 


The galactic potential <I> is appropriate for a stellar density distribution p*(r) and a dark matter halo. 
Except for the magnetic force term in the equation of motion (3.2) and heating due to reconnection of 
magnetic fields, the last term in equation (3.3), these equations are identical to those used by Brighenti and 
Mathews (1996). 


In the absence of plasma resistivity and field sources, the magnetic field evolves according to the usual 
frozen-in condition: 

— = Vx(vxB). (3.4) 

However, for applications to the interstellar medium of elliptical galaxies, it is necessary to include source 
terms due to fields ejected by evolving stars, turbulent growth, and the loss of field due to reconnection. 
These modifications will be developed in the following sections. 

Although the non-magnetic source terms in equations (3.1) - (3.3) have been described in detail by 
Brighenti and Mathews (1996), for completeness we shall briefly review them again here. New gas is 
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introduced into the interstellar medium by mass ejection from normally evolving stars formed in a single 
burst at time t = 0. The specific rate of mass loss is proportional to the local stellar density (gm cm“^ 
s“^) where 

(3.5) 

(Mathews 1989). Here = 15 Gyr represents the present time and a(t„) = 5.4 x 10“^° s“^, valid for a 
wide variety of initial mass functions (Mathews 1989). On average this new gas appears in the galaxy at 
zero velocity relative to the net local stellar velocity v* which is non-zero only for rotating galaxies. The 
negative source term in equation (3.2) represents the momentum drag on the local gas flow v as new gas 
created with velocity v* is accelerated to the local flow velocity v. 

Equation (3.3) for the conservation of specific thermal energy e = iksT/2^mp (nip is the proton 
mass) contains a negative term representing optically thin radiative losses by thermal emission as computed 
by Raymond, Cox & Smith (1976): nen//L(T)/p = fpL{T)/rn^ ergs s“^ gni“^ where / is related to the 
mean molecular weight by / = (2 -|- /r)(4 — 3^)/25^^; we shall assume p = 0.5. Additional source terms in 
equation (3.3) describe the heating of the gas by the thermal and kinetic energy of stellar ejecta, including 
supernova explosions. After new gas enters the flow with mean velocity v* it is accelerated to the local gas 
velocity, dissipating energy at a rate Q!*/3 *|v — v*p/2 per unit volume. The dissipational heating resulting 
from random stellar motions and supernova explosions is represented by a mass-weighted characteristic 
temperature for the new gas, 

To = (a*T* -I- as„Ts„)/a (3.6) 

where « = «*-!- Ugn ~ cc* and Eq = 2>kBTo/2pmp is the corresponding specific thermal energy. The term 
P!p in the source term of equation (3.3) represents the work done by the stellar ejecta on the local ambient 
gas as it expands from an initially dense cloud toward pressure equilibrium with the interstellar gas. 

The solar or subsolar iron abundance observed in the interstellar medium of ellipticals indicates that 
the (Type la) supernova rate is low and that it has been low in recent Gyrs (Loewenstein and Mathews 
1991; Loewenstein et al. 1994). Observations of the current supernova rate in elliptical galaxies also suggest 
rather low values, ~ 0.1 SNu (van den Bergh, McClure, & Evans 1987; Cappelaro et al. 1993) where 1 SNu 
corresponds to one supernova per 10^*^ Lb every 100 years. The past rate of supernova explosions is even 
more uncertain, but it is generally assumed that the rate was higher when the galaxy was younger. We 
shall assume a power law time dependence of the form 

asuTsu = 2.13 X 10-® SNu Ks”! (3.7) 

where Lb is in units of the solar value Lqb = 4.98 x 10^^ ergs s“^ and the reduced Hubble constant is 
assumed to he h = 77/100 = 0.75. The SNu parameter defined by equation (3.7) cannot be compared 
directly with observed values since other uncertain parameters (such as Egn) also appear in this same 
equation. Recent calculations of spherical cooling flows (Mathews 1997) have shown that global thermal 
instabilities or “galactic drips” can occur in isolated ellipticals if the supernova rate parameter SNu ^ 0.05. 
Therefore we shall assume SNu = 0.066 and Esn = 1 here although the iron production may still be 
somewhat larger than observed (Mathews 1997). The time variation for asnTsn oc is similar to that 
adopted by Ciotti et al. (1991) but less steep. 
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4. MAGNETIC COOLING FLOWS IN SPHERICAL GEOMETRY 

4.1. Galactic Model 


We choose a King model for the stellar density 

p*ir) = po*[l + 

where Tc* is the stellar core radius. Most of the gravitational potential arises from a dark halo having an 
approximately isothermal distribution: 


Ph{r) = poh[l + [r/rchf] ^ 

Parameters describing the spherical galaxy model listed in Table 1 have been chosen so that the galaxy lies 
on the fundamental plane (see Tsai & Mathews 1996a and Brighenti & Mathews 1996 for details). The 
stellar temperature T^:(r) = {pmp/3k)a^ can be found by solving the equation of stellar hydrodynamics 
in spherical symmetry (Binney & Tremaine 1987; Mathews 1988). For this purpose we assume a velocity 
ellipsoid of the form 

(V^r) - {VD 

(v?) 

where and (u^) are the transverse and radial stellar velocity dispersions respectively. We choose q = 2, 
but the evolution of the ISM is not sensitive to the value of q over a wide range 1 ^ q ^3. 



4.2. Spherical Magnetic Cooling Flows 


In spherical symmetry, the magnetic field can have radial and tangential components, B = {Br,Bt) 
and equation (3.4) becomes 

d log B^ 4u 

dt r 

and 


(4.1) 


d log Bf 
dt 


r 


d\ogu 

d\ogr 


(4.2) 


where d/dt is the Lagrangian or comoving derivative and u(r) < 0 is the radial velocity of the cooling flow. 
The (source-free) equation of continuity has a similar form: 


dlogp d\ogu 

dt r \ i91ogr^ 

Dividing equations (4.1) and (4.2) by (4.3), we find 

dlogi?^ _ 4 _ 

dlogp ( 2 -^ 7 )“"^’' 


(4.3) 


(4.4) 


d\ogB^ _ 2(1-by) ^ 
dlogp (2-by) 


and 


(4.5) 



where 


(4.6) 


_ 91ogM 
i91ogr 

Therefore each component of the field evolves with the local gas density as a power law having an exponent 
that depends on the structure the local velocity field as characterized by the exponent 7 . In general 7 can 
be a function of galactic radius. 


For example in a homologous, Hubble-type flow, u cc r/t and p is spatially uniform. In this case 7 = 1 , 
7r = 7i = 4/3 so the total magnetic field 


=Bl + Bf 


(4.7) 


evolves according to 


dB"^ _ 4 ^^ 
dt 3 p dt 


(4.8) 


In typical galactic cooling flows the gas velocity slowly accelerates toward the galactic center, |m| cx l/r^ 
where 0 < p < 1. For this range of velocity fields the radial component B,. grows faster than Bt toward the 
(high density) center of the cooling flow (e.g. i?/ oc Bf ~ constant for p = 0; Bf gc p^, Bf (X p ior p = 1). 
As a result of this unequal growth, an initial field random on small scales in non-rotating ellipticals will 
develop radial fields having large scale coherence as the flow approaches the galactic core. Differential shear 
in slowly rotating ellipticals is expected to convert disordered fields into coherently toroidal configurations 
(Brighenti & Mathews 1996). In either case, however, reconnection may be expected to occur in regions 
of random field; when such disordered fields are stretched by the global flow, reconnection could still be 
important since field lines in adjacent flux tubes will tend to be antiparallel. 


The fraction of the field contained in random and coherent (but antiparallel) components in galactic 
cooling flows can be estimated by comparing the bulk flow velocity u and the mean turbulent velocity Vt', 
both velocities are subsonic. Whenever B^ is less than the equipartion held Bf and u ^ Vt we expect 
that the field geometry will be continuously randomized by local turbulence; should u > Vt the field will 
be stretched with the fluid and develop large scale coherence. The ratio of characteristic timescales for 
turbulence and bulk flow tturb/tfiow = ufvt can be used to define a parameter 


6{r) = 1 — expiy—u/vt) 


that increases with the development of large scale fields. Here 5 ^ 0 or ^ 1 in strong or weak turbulence 
respectively. In general we may imagine that the individual field components have a (source-free) evolution 
described by 


dBf 

dt 




and 


dt 


= [(57t-b(l-(5)^]Hj 


2 

p dt 


When continuous randomization obtains, <5 « 0 implies B^ oc p"*/^ and equation (4.8) holds; in the following 
discussion we shall assume this limit, i.e. the field is randomized on timescales short compared to the (quite 
long) flow time in the interstellar medium. The validity of this assumption can be estimated by evaluating 
S{r) from the computed cooling flow solutions. 



- 9 - 


4.3. Seed Fields, Turbulent Growth and Reconnection 


While the flux-freezing condition (4.8) describes the evolution of a completely disordered field, 
additional source and sink terms must be included to account for the origin, turbulent growth and loss of 
the field: 


dB^ 

dt 


3 p dt Tturb 


bD- 




+ 


op* 


Bi 


+ CC7m-87rP(5(t -tjn)- 


(4.9) 


The second term on the right represents the growth of the field due to turbulent amplification; the unit 
step function [9{x) = 1 a; < 0; 9{x) = 0 a; > 0] ensures turbulent growth only when B^ < B^ = AnpVf. This 
term deals in a schematic fashion with a subject of considerable complexity. Normally the turbulent growth 
of the held would be represented using the (linear) mean held dynamo equation (see Beck et al. 1996 and 
references therein), however in view of the difficulties described earlier in generating equipartition with 
the mean field equation, we prefer the simplified approach in equation (4.9). The third term on the right 
accounts for the loss of field by reconnection. The last two terms on the right represent original seed fields 
from stellar mass loss and galactic events at time tm which create non-stellar turbulence in the interstellar 
gas. These terms are discussed in more detail in the following. 


4-3.1. Seed Fields 


Aside from the possibility of primordial fields which are very small in any case (Beck et al. 1996) the 
most credible source of magnetic field in the interstellar medium are the fields created by dynamos in stellar 
interiors which are expelled into the ISM when red giant stars lose their envelopes. The relevant magnetic 
fields are not those observed in stellar photospheres or chromospheres, but the mean fields that are 
likely to exist throughout (convective) red giant envelopes before ejection from the star. Guided by the 
magnetic flux observed in white dwarfs. Moss and Shukurov (1996) adopted 6* = B^jp^J^ « 2 x 10® in cgs 
units as a typical value, corresponding to B^ « 10^ G and p* w 5 x 10® gm cm“®; we adopt this same value 
here. After the stellar envelope has expanded to the local density of interstellar gas n = p/rup conserving 
magnetic flux, the seed field will be very small: Bg = « 2 x 10“^®(n/10“®)^/® G. 

An additional source of seed fields may arise from strong turbulent mixing following a significant merger 
of another galaxy with the elliptical at time tm or after some other global environmental disturbance. 
Following such an event, we assume that the entire ISM becomes strongly turbulent, amplifying any seed 
fields until the magnetic energy density is some fraction Wm of the thermal energy density in the gas. Such 
a hypothetical turbulent event is represented schematically in equation (4.9) by a magnetic energy density 
proportional to the local gas pressure, i?^/87r w n7m(3P/2) imposed suddenly at time tm- 


4-3.2. Turbulent Growth of Field 


The growth of magnetic flux by plasma turbulence in the interstellar medium of elliptical galaxies has 
been discussed in detail by Moss and Shukurov (1996); we shall adopt many of their suggestions here. The 
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two principal internal sources of interstellar turbulence are the mixing of gas ejected from orbiting stars and 
the creation and buoyant evolution of (Type la) supernova remnants (Mathews 1990). We assume that the 
random vortical motion associated with these sources amplifies the interstellar magnetic field. The energy 
density in the field is expected to grow until it becomes comparable with that in the turbulence and then 
saturates at that value. The equipartition field is defined by 

1 „ 3 3 , , , 

= ^^PCs- (4.10) 


The turbulent energy density is written as a fraction tz7(r) of the local thermal energy density (c^ is the 
isothermal sound speed and Vt is an rms speed characterizing local turbulent motion). Normally, we 
expect tz7 <C 1 and therefore ^ 3Pj2. Ideally, a separate dynamical equation should be introduced 

to describe the creation, transfer and evolution of turbulent energy density. However, we avoid these 
complications here by assuming that the turbulent energy density is generally small {zu < 1) and that the 
principal catalytic role of turbulent energy is to convert sources of interstellar energy into magnetic energy 
density. In addition we shall assume for simplicity that the turbulence generated by stellar ejecta and 
supernovae is homogeneous and isotropic; this (somewhat optimistic) assumption is discussed by Moss and 
Shukurov (1996). 


The turbulent field is produced by a superposition of kinetic energy created by stellar ejecta and 
supernova remnants. The stellar contribution to the local turbulent energy density is given by the product 
of the rate of energy production by mass ejection from evolving (red giant) stars Q;*/9*(3A:T*/2/imp) (ergs 
cm“^ s“^), the fraction Ct of this energy that goes into turbulence, and the characteristic eddy turnover 
time r = itjvt: 



a^,p^,et 


SkT^, £* 

2p,mp 




(4.11) 


The characteristic rms turbulent velocity is = (3cc7*)^/^Cs although the full range of turbulent velocities 
must extend up to the local stellar velocities.. The typical eddy size should be comparable to the 
dimension of a stellar envelope (or planetary nebula) after it has expanded to pressure equilibrium (at 
T « lO^K) with the local interstellar gas. For this purpose (following Moss & Shukurov 1996) we use the 
estimate of this length scale determined by Mathews (1990), « 4.7 x 10®/?“^/^ cm, where p is the local 

interstellar density. We assume that the turbulent mixing time is the eddy turnover time of the largest 
eddies, this is conservative since the turnover rate is faster for smaller eddies. Equation (4.11) can 

be solved at each galactic radius for the stellar equipartition parameter 


OT*(r) 


1 

3c2 



2/3 


(4.12) 


where cl = kT^/pmp. For the elliptical described in Table 1 T* w 2.95 x 10® K (for 10® ® ^ r^pc ^ 10^-^®) 
and p* w 4.35 x 10“^^r)rp® gm cm“®. Field-free cooling flows in this galaxy result in gas density and 
temperature distributions given by p « 1.05 x 10' "^^Tkpf gm cm ® and T « 8.91 x 10®rj,j)(,'^® K at time 
t = 15 Gyr. These relations result in a stellar equipartition parameter n7* w 1.07 x 10“"‘r^p®^'^®e(^^(t/t„)“® ®’^ 
valid for 10® ® < rkpc ^ 10^-^®. The corresponding turbulent velocity induced by stellar ejecta, 

« 6.9e(^^r^p®g'^®(t/t„)“® "‘® km s“^, is indeed very subsonic (c* « 385r^p®g'^® km s“^). The turnover time 
for the largest eddies is r* = ~ 4.6 x lO'^e)”^^^?’®^® (t/t„)“®-^® years. 

In a similar manner Type la supernovae feed energy into the interstellar medium at a rate 
Q;*p*(as„Ts„/a*)(3A:/2/tTOp) (ergs cm“® cm“^) and a fraction et of this is assumed to go into turbulent 
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motions. The eddy mixing time is tsnivsnt where Vsnt = (3i37sn)^^^Cs. For supernova-driven turbulence the 
(largest) eddy dimensions are comparable to the size of the hot bubble formed in the interstellar medium 
following each supernova event, ^sn = 3.6 x 10^^/?“^/^ cm, based on the blast wave calculations of Mathews 
(1990). Combining these expressions, the equipartition parameter for supernova induced turbulence is 


Wsnir) 


1 { p* ik asnTsnisn \^^^ 

3 V * P MTOp c 3 J 


(4.13) 


For the bright elliptical shown in Table 1 we Hnd Wgnir) ~ 1.11 x 10“^ (SNu/0.01)^/^(t/t„)“^/^r^p°’^^ 

(valid for 10°'^ ^ Vkpc ^ lO^ "^^), which is somewhat larger than the stellar contribution to turbulence 
for a given e*. The corresponding rms turbulent velocity and turbulent mixing time are Vsnt ~ 
22.2ei/^rfcp°-2°(t/t„)-i/6(SNu/0.01)i/3 km s'^ and n = e*/v*t « 1.1 x 10V^^^^^°(V^n)^/®(SNu/0.01)-i/3 
years. 

The combined turbulent mixing time required in equation (4.9) is found from 


1 _ ^ '^snt 

'^turb ^sn 


(4.14) 


and the total equipartition field is given by equation (4.10) with va(r) = zu^.ir) + vasn{i~)- Obviously the 
introduction of turbulent energy by supernovae occurs at the location of these events, not in a spatially 
smooth manner as represented by the turbulent source term in equation (4.9). Nevertheless we shall assume 
that the turbulence is driven in a smooth, homogeneous manner. The applicability of this assumption is 
discussed by Moss & Shukurov (1996). 

In addition to using stellar processes to produce vortical turbulence, Moss and Shukurov (1996) also 
invoke “acoustic turbulence,” a general field of compressional (sound) waves or noise associated with stellar 
mass loss and (particularly) supernova explosions. Although the coupling of this acoustical energy density 
to the usual vortical turbulence is non-linear and inefficient, Moss and Shukurov claim that acoustical 
wave energy contributes substantially to the turbulent growth of the magnetic held and that these waves 
ultimately steepen into shocks. However, in the absence of reasonably strong magnetic Helds, a condition 
that may obtain at early times, viscosity and thermal conductivity in the hot interstellar medium are 
expected to damp acoustic waves within a wavelength of their source. Compression waves should damp 
before they steepen into shocks. For example, damping due to plasma viscosity p exponentially reduces the 
flux of plane parallel sound waves of frequency uj = ^ixcsjl propagating in the x-direction by where 




= 6 X IQi^ 


^pc 


( _^) 

V 10“^cm“3 J 



cm 


(e.g. Zel’dovich and Raizer 1966). Damping by thermal conductivity k is comparable. Evidently acoustic 
waves generated by either stellar mass loss or supernova events are strongly damped and thermalized within 
a single wavelength of their point of origin in the interstellar medium. Both p and k. may be reduced by the 
presence of a magnetic Held, but a reduction factor of ~ 10^ would be required for acoustic damping to be 
negligible, requiring the presence of signiHcant magnetic Helds at early times. In view of this complication 
we do not consider turbulent Held generation by the acoustical wave Held here. 
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4-3.3. Field Reconnection 


The expected dominance of random magnetic fields in elliptical galaxies implies that the occurrence of 
counter-directed nearby fields should be common. In such regions thin boundary layers are created where 
the finite resistivity converts field energy to other forms of energy: thermal, bulk kinetic energy, accelerated 
particles, etc. (Soker & Sarazin 1990; Lesch & Bender 1990; Jafelice & Friaca 1996). The physics of field 
reconnection in the hot ISM of ellipticals galaxies is discussed in some detail by Jafelice and Friaca where 
many earlier references are cited. From dimensional and physical arguments, the flow into the reconnection 
regions is expected to occur at some fraction of the local Alfven velocity Va- We shall follow the simple 
representation used by these authors for the field reconnection time scale in equation (4.9): 


'^rec 


£rec _ ^rec(47rp)^/^ 
Va B 


(4.15) 


Here ^rec is a single parameter that describes the reconnection efficiency. In addition we shall assume with 
these authors that all of the field energy lost by reconnection is entirely converted to thermal energy; so a 
term /SiTTrec must appear in the thermal energy equation (3.3). The principal parameters related to the 
turbulent dynamo are et and Irec both of which are very uncertain. 


4.4. Spherical Cooling Flow Equations 


Summarizing the previous results, the equations describing the time-dependent evolution of spherical 
cooling flows containing a self-generated, disordered magnetic field are: 


du 


dp 

m 

dP 

dr 


1 d{r‘^pu) 

r 2 Qj, 

d(rB^) 


1 


127rr dr 


= ap*, 
9$ 

dr 


api,u, 


de 


P dp 
p dt 


P^fL{T) 


Otp^: 


P 1 2' 

Sot -£ + -zU 

P 2 _ 


B^ 


STTTr, 


The source of thermal energy for the gas is determined by 


(4.16) 

(4.17) 

(4.18) 


Sot = eo[l - et0(H2 - Bl)] 


which allows a fraction (1 — et) of the energy from stellar processes to be thermalized when B < Be- When 
B > Be all the stellar energy is assumed to heat the gas either by direct heating or by eventual dissipation 
of the turbulent component. The magnetic field evolves according to equation (3.9) which we write as 


dB^ _4B^dp 
dt 3 p dt 


dB^ 

dt 


I 


B^ 

Vrec 


ap* 

p 


^ j p'^^^ + Wm'^SzP5{t - tm). 


(4.19) 


We consider two simple representations {i = 1,2) for the turbulent growth of the magnetic field: 

dt J ^ Tturh 


-0{B^ - Bl) Case I 


(4.20a) 
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and 

= 8TTap^:eoetO{B‘^ — b1). Case II (4.206) 

2 

For Case I the magnetic field can increase only on the time scale Tturb that characterizes turbulent growth. 
The step function 9{B'^ — B^) guarantees that the field cannot grow beyond equipartition. In Case I some 
part of the power density ap^,eteo is allocated to increase the field strength and the rest is assumed to 
maintain the local turbulent energy density. 

Alternatively, the Case II option allows stellar and supernova energy to convert directly to either 
magnetic or thermal energy density, depending on the value of 6{B'^ — i?g). The turbulent energy density, 
which still can be regarded as an intermediate stage, is not explicitly considered. The field grows at a rate 
87 rap*et£o until it reaches local equipartition then growth suddenly saturates. When equations (4.18) and 
(4.19) are added in the Case II approximation, the total energy is explicitly conserved either as thermal or 
magnetic energy. In the following we consider solutions to both sets of equations. Cases I and II. 



5. Cooling Flow Models 


We have computed evolving magnetic cooling flows in the galaxy described in Table 1 using a 
Lagrangian code. The primary parameters that define each flow are SNu, e* and irec- Our cooling flow 
solutions are not sensitive to the value of the seed field parameter 6 * = B^/pt^^ since for all cases considered 
the exponential growth of field in the turbulent dynamo is rapid compared to the local gas flow time (~ r/u) 
in the ISM. The parameters zum and tm (see eqn. 4.19) can be used to simulate the global turbulence 
expected following a galactic merger or some other short-lived source of turbulence. Each calculation is 
begun at time t = 1 Gyr when it is assumed that the galactic winds (driven by Type II supernovae) have 
reversed to become cooling flows (David et al. 1990; 1991). Interstellar turbulence is turned on at a later 
time (6 = 2 Gyr) to avoid numerical interference that may occur if the gas and the magnetic field are 
initiated at the same time. The calculations are stopped at 15 Gyr which we regard as the present time. 
Our elliptical galaxy is assumed to be perfectly isolated apart from the environmental influences represented 
by the parameters Wm and tm- 

In Figure 1 we compare several cooling flow models at time 6 = 15 Gyr. The non-magnetic flow shown 
in Figure la-c resembles spherical cooling flows that have been discussed in detail elsewhere (e.g. Brighenti 
& Mathews 1996; Mathews 1997). The parameters used for this non-magnetic cooling flow are: e* = 1 
and SNu = 0.066. Some numerical irregularity is visible in Figures la and lb for logr ^ 2.5 because of 
the smaller number of computational zones in the galactic core. Most of this numerical noise arises when 
central zones are removed as they cool below T = 10'^ K. We recognize that the x-ray surface brightness 
Sa;(r) corresponding to our model in Figure la-c decreases somewhat more steeply with galactic radius 
than that typically observed (Trinchieri, Fabbiano, & Ganizares 1986). The traditional correction for this 
discrepancy has been to remove gas from the flow, assuming that stars or some other invisible objects are 
formed from the mass that “drops out” (e.g. Stewart et al. 1987; White & Sarazin 1987; Thomas et al. 
1987). Our models would adjust in a similar manner if mass were removed. We have not included drop out 
in our calculations for several reasons: (i) we wish to illustrate the possible dynamic influence of magnetic 
forces without additional uncertain complications, (ii) in those (perhaps unlikely) situations when the 
magnetic field is strong, its influence on the solution (and '^^{r)) is very similar to mass drop-out, and (iii) 
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thermal instabilities, often invoked to justify mass drop out, are generally arrested at a very early stage in 
the presence of modest magnetic fields. 

The corresponding Case I magnetic cooling flow evolution with parameters Ct = 0.5 SNu = 0.066, 
^rec = 3 X 10^° cm, and zum = 0 is illustrated by the solid lines in Figure Id-f. Our choice of a large ^rec 
effectively shuts down the reconnection process, maximizing the computed field strength. Even though half 
of the energy from stellar sources is available to drive the turbulent dynamo, the gas density, temperature 
and pressure are all essentially unchanged from the field-free cooling flow (Figure la-c). The magnetic 
pressure (long dashed line in Fig. If) is never more than a few percent of the gas pressure although 

it is somewhat larger at small galactic radii. But the field strength indicated in Figure If (long-dashed 
curve) is not small, reaching 1 /iG at r = 7 kpc and 10 /rG at 1 kpc. Of course these fields are upper limits 
because we assume that the dynamo is efficient and our choice of large free effectively deactivates field 
reconnection. We find that the results in Figure If are essentially unchanged as long as Gee ^ 10^® cm. 
The magnetic field throughout the galaxy slowly decreases with free; in Figure If we show that the field 
is reduced only by a factor of ~ 2.5 (short-dashed line) when free = 10^^ cm. Unfortunately the physical 
significance of this is unlear. The reconnection process, if it is important, is expected to occur throughout 
the flow in very thin boundary layers separating regions of adjacent antiparallel fields. Such a physically 
complex process cannot be adequately modeled with the single parameter free- Also illustrated in Figure If 
(dotted line) is the magnetic pressure at f = 15 Gyr that results when et is increased to 0.8 with the other 
parameters unchanged (SNu = 0.066, free = 3 x 10^° cm, and Wm = 0). The field is slightly larger in this 
case because of the somewhat greater turbulent energy that existed early in the cooling flow evolution. 

We have also computed cooling flow solutions with the same parameters but using the Gase II 
equations. We find that the results are essentially identical to those of the Gase I solutions. As explained 
below, the insensitivity to the specific representation (Equation 4.20a or b) for the turbulent growth of the 
field arises because for either Gase the field grows to equipartition early in the calculation in a time short 
compared to the local flow time (~ r/u). 

Of particular interest for understanding the amplification of disordered fields in galactic cooling flows 
are the radial dependences of the thermal, magnetic and turbulent energy densities. Using the approximate 
power law variations discussed in §4.3.2, the thermal energy density should vary as Ctherm oc pT oc 
Eor turbulence driven by either stars or Type la supernovae, the turbulent energy density varies as 
^turb = pVt/‘^ c)c ptuCg oc and the equipartition field must vary in the same manner, oc . 

However, the field compression that accompanies the cooling flow causes a slightly steeper variation, 
emag OC OC 

Two important conclusions follow from these spatial dependencies: (1) the increase in the magnetic 
field strength toward the galactic center can be dominated by compression in the cooling flow, not local 
turbulence, and, if so, (2) the magnitude of the interstellar field depends entirely on the most distant radius 
in the galactic ISM where the turbulent dynamo process amplified stellar seed fields at some previous time. 

In Figure 2a we plot the spatial variation of the thermal, magnetic and turbulent energy densities 
for the standard magnetic cooling flow solution (shown in Figure Id-f with solid lines) at times t = 12 
and 15 Gyr; these contours are essentially identical for Gase I and Gase II equations. The slopes of the 
various energy densities (in the linear parts of the log-log contours) are in excellent agreement with the 
approximate variations estimated above. At either 12 or 15 Gyrs tmag is almost 10 times larger than eturb 
throughout most of the galactic volume. As expected from the previous estimate, the local turbulent energy 
density plays no explicit role in field amplification since > B^. This excess of the field above the local 
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equipartition value occurs because the field at t = 15 Gyr is a relic of a historical field which developed 
when the turbulence was larger. At an early time in the evolution of this cooling flow equipartition 
(cmag ~ ^turb) was reached in time ~ T^urb throughout the galaxy, but as the calculation proceeded the 
turbulent energy decreased {tturb oc oc as expected from the scaling relations in §4.3.2 when 

vlsn ~ ’v't*)- The super-equipartition held is maintained and amplified by compression in the cooling 
flow. Since eturb/^mag increases with radius, the current galactic field is most sensitive to past turbulent 
conditions in the outermost parts of the galactic flow. 

In this regard we note the small decrease in Cmag in Figure 2a that occurs near the outer boundary rt 
of the galaxy between 12 and 15 Gyrs. In this part of the flow the gas density drops slightly during this 
time interval due to the (losing) competition between the decreasing rate that gas is introduced by evolving 
stars and the depletion of gas by global inflow toward the galactic core. As p decreases, the magnetic field 
must also decrease oc Some additional loss of local field near may occur because the local 

field is being diluted by new gas that only contains the very small seed field yet turbulent amplification 
is not possible until drops below the equipartition value. Figure 2b shows how the frozen-in field has 
been advected inward toward the galactic center, creating and maintaining the expected slope oc 
throughout the interstellar medium. This also illustrates how the field throughout the X-ray bright central 
parts of the ISM depends critically on the historical turbulent growth in very distant regions. In view of 
the relationship between the global field and past turbulence near the outer edge of the galaxy, it might 
be thought that the growth of the field would be sensitive to the outer boundary conditions in the flow, 
but this is not the case. For example we repeated the standard calculation shown in Figure Id-f (solid 
lines) with an open boundary condition at r = which allows gas to flow out beyond the galaxy into the 
surrounding vacuum. After t = 15 Gyrs the flow variables and the field strength distribution within the 
galaxy were essentially identical to those shown in Figure Id-f. Evidently the gas density and pressure near 
r < rt is maintained by local stellar mass loss, not by the small amount of gas that flows out of the galaxy 
when an open boundary is allowed. 

Finally, we have performed some calculations with non-zero Wm to simulate an event, such as a galactic 
merger, that introduces global turbulence throughout the ISM at some time tm in the past. In Figure Ig-i 
we illustrate the magnetic cooling flow at time t = 15 Gyr following a disturbance in the the interstellar gas 
that occurred at tm = H Gyr when the turbulent energy density was raised to one fourth of the thermal 
energy density {wm = 0.25) throughout the galaxy. As a result of this event, the magnetic energy density 
at 15 Gyr (Fig. li) is ~ 100 times larger than it would have been otherwise (Fig. If). This illustrates 
that present day magnetic fields in elliptical galaxies can be significantly increased by past events that 
stirred up the outer parts of the interstellar gas and which expended energies that are small compared to 
the total binding energy of the interstellar gas. For the cooling flow in Figure Ig-i the magnetic energy 
density zUm3P/2 = 3P/8 is imposed impulsively at time tm and subsequent growth of B'^IStt is due to 
compression in the global flow. By time t = 15 Gyr Figure li shows that magnetic pressure dominates the 
flow dynamics within about 5 kpc of the galactic center while the gas density remains monotonic, implying 
that the magnetic core does not become buoyant. However, the magnetic stresses are sufficient to slow 
down the local flow in r 5 kpc allowing an accumulation of gas near this radius where some radiative 
cooling and mass deposition occurs. This is reminisent of the “galactic drip” phenomenon in which gas in 
non-magnetic cooling flows can cool far from the galactic center (Mathews 1997), but the increase in gas 
density and thermal instability shown in Figure li has a purely magnetic origin. When Wm is lowered to 
0.15, B^jSir w P occurs at a smaller radius, r « 1 kpc, and radiative cooling beyond the galactic core no 
longer occurs. 
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Some caution must be exercised in interpreting the non-stellar turbulent growth imposed at tm in 
terms of a galactic merger. For galactic mergers in which the galactic mass is increased by a relatively small 
factor, ~ 10 percent, Mathews and Brighenti (1997) have shown that enough energy can be released by 
Type II supernovae from newly created stars to remove most or all of the preexisting galactic interstellar 
gas. Such an energy release obviously depends on the assumed initial mass function for the newly formed 
stars. We do not consider these complications here. 


6. DISCUSSION AND CONCLUSIONS 


In this study of the evolution of interstellar magnetic fields in elliptical galaxies we have shown that 
the observed fields can be understood as a natural result of internal evolutionary processes: stellar seed 
fields are first amplified by a turbulent dynamo then by compression in the galactic cooling flow. We have 
assumed that stellar dynamos occur in virtually all stars and that these stars supply tiny seed fields to the 
ISM as they undergo mass loss during normal red giant evolution. Our optimism concerning the efficacy 
of the turbulent dynamo mechanism and the relative inefficacy of field reconnection have been essential in 
obtaining central field strengths comparable to those observed, ~ I — 10 fiG. 

Our models extend a modified version of the turbulent dynamo theory proposed by Moss and Shukurov 
(1996) to include the important effects of compressional field amplification in the global galactic cooling 
flow and the evolutionary relationship between currently observed fields and those created by galactic 
turbulence many Gyrs ago. As the field advects from the outer regions of the galaxy to the galactic core 
the field compresses by > 10^ and the magnetic energy density increases by > 10®. In our models the 
magnetic energy density exceeds the local turbulent energy density {B^ > Bf) throughout most of the 
interstellar gas so local turbulence plays little or no role in field amplification. When B^ > B^ the field may 
quench some of the turbulent activity although we have not allowed for this feedback in our models; it is 
also unclear whether the random character of the field (on scales of t'* or £sn) can be fully maintained in 
this limit. For most of the magnetic cooling flows described here we have assumed that half of the stellar 
energy goes (at least initially) into turbulence (e* = 0.5). Although the ratio et of turbulent energy to 
immediately thermalized energy resulting from stellar mass loss processes is uncertain, our results suggest 
that the resulting field strength does not depend strongly on et- 

How sensitive are our results, particularly the condition B^ > B^, to the specific galaxy model we 
have adopted, a simple King stellar distribution with an isothermal dark halo? To address this question 
it is useful to express the radial variation of the magnetic and turbulent energy densities in terms of the 
gas and stellar densities. For either stellar or supernova turbulence we find Cturb = pVf/2 (X and 

Cmag oc (X Therefore the ratio of these energy densities Cmagliturb oc Clearly the 

radial variation of the gas density depends on both the stellar and dark matter density profiles, so the 
interpretation of this last relation is complex. To make further progress, suppose that we adopt the point 
of view that the X-ray and optical surface brightness profiles are the same, i.e. E* oc Canizares, 
Fabbiano, & Trinchieri (1987) found this to be the case for three bright ellipticals in Virgo, but more recent 
ROSAT observations show some deviations from exact proportionality. However, if S* oc then p^ oc p* 
and emagl^turb OC SO that €mag increases very slowly relative to Cturb with increasing galactic radius, 

opposite to the variation shown in Figure 2a. We may expect, therefore, that the ratio Cmag/^turb could 
either slowly increase or decrease with galactic radius, depending on the details of the galactic model. 
However, regardless of these details it is clear that the interstellar magnetic field can increase toward the 
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galactic center no slower than oc as shown in Figure 2b. Finally we note that in our models the 
^magl^turb increases with the cooling inflow most strongly in the outermost parts of the galaxy where the 
power law scaling breaks down and where the observational uncertainty is greatest. 

If interstellar fields in elliptical galaxies have their origins in the seed fields of stellar ejecta, the 
physical scale d of regions of Faraday coherence may be smaller than previously supposed (Burn 1966; 
Carrington & Conway 1991). Regions in the hot interstellar gas having masses comparable to those of 
stellar envelopes, 0.4 Mq, have diameters of only d « 30(n/10 x 10“^)“^/^ pc. However, Faraday-coherent 
regions of parallel field lines will be even smaller because of the violent instabilities that accompany stellar 
mass loss. Since the observed Faraday depolarization parameter is proportional to (ni?|| the field 

implied by any observation will be larger as d is reduced. This difficulty is not alleviated by large scale 
fields that can develop in (non-Hubble) spherical or differentially rotating non-spherical cooling flows, since 
such an “aligned” field will still be counter-directed on small mass scales. The scale of field coherence d 
can increase if reconnection is efficient, but this would also decrease the overall magnetic energy density. 
Clearly, more attention must be given to the likelihood of field reconnection in cooling flow environments. 

Our main conclusions are: 

(1) Self-generated interstellar magnetic fields in elliptical galaxies can be comparable to those suggested in 
Faraday depolarization studies provided (i) the turbulent dynamo process can amplify the mean field to 
equipartition with the turbulent energy density and (ii) field reconnection is not dominant throughout most 
of the ISM. 

(2) Using our galactic model, a significant additional amplification of the overall field occurs by compression 
in the galactic cooling flow; the magnetic energy density is increased by factors of > 10®. 

(3) Throughout most of the galactic volume the field energy density can exceed local equipartition with the 
turbulent energy density. 

(4) The exact value of the field in ellipticals is sensitive to interstellar turbulence at large distances from 
the galactic center and at distant times in the past. For this reason we may expect some variability in field 
strengths among elliptical galaxies, depending on the environmental history of the outer galaxy. 

(5) In our spherical cooling flow models the magnetic energy density or pressure in self-generated cooling 
flows {zum = 0) never dominates the galactic gas pressure (nor becomes buoyant) so we do not expect the 
X-ray images to be strongly influenced by the field. In more realistic rotating cooling flows, however, the 
magnetic pressure resulting from differential rotation may decrease the X-ray brightness near the cooling 
flow disk (Brighenti & Mathews 1996). 

(6) If the interaction of the elliptical with ambient gas or nearby galaxies stimulates strong turbulence in 
the galactic ISM, the resulting field strength may be greatly enhanced many Gys afterward. 


We are pleased to acknowledge comments from Niel Turner and support from a UCSC Faculty Research 
Grant and NASA grant NAG 5-3060. 
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Table 1. GALACTIC PARAMETERS FOR ELLIPTICAL GALAXY 


Parameter 

Value 

Tc* 

311.59 pc 

^ a 
' e 

5.088 kpc 

Teh 

6.22 kpc 

rt 

113.1 kpc 

P^o 

1.438 X 10“^° gm cm“^ 

Pho 

5.46 X 10“^'^ gm cm“^ 

Mift 

4.52 X 10“ Mq 

Mht 

9 M^t 

Lb 

4.95 X 1O1°Lb0 

Mtft/L b 

9.14 

CT*'’ 

351 km s”'^ 

log 

40.0 


^'Effective radius. 

'^Characteristic velocity dispersion in stellar core, cr* = 
(47rGp*or2*/9)^/^. 

'^For 0.5 - 4.5 keV based on correlation of Donnelly, Faber 
& O’Connell (1990); more recent values are higher (Kim & 
Fabbiano 1996). 







Figure Captions 


Figure 1: Three vertical arrays of figures show cooling flow conditions in the interstellar gas after evolving 
to t = 15 Gyr. Solid lines in each plot refer to the gas density, temperature or pressure corresponding to 
these parameters: (a-c), nonmagnetic cooling flow with e* = 0.5 and SNu = 0.066; (d-f), magnetic cooling 
flow with €t = 0.5, SNu = 0.066, £rec = 3 x 10^° and zum = 0; (g-i), magnetic cooling flow with e* = 0.5, 
SNu = 0.066, irec = 3 X lO^o, 'cc'm = 0.25 and tm = 11 Gyr. Long-dashed lines in the density plots show 
the variation of the (unnormalized) stellar density. Long-dashed lines in the temperature plots show the 
variation of stellar temperature. Long-dashed lines in the pressure plots show the magnetic pressure. Plot 
(f) also shows the magnetic pressure with tt = 0.8 (dotted line) and with some field reconnection Gee = 10^^ 
(short-dashed line). See text for further details. 

Figure 2: (a) Variation of energy densities with galactic radius at t = 12 Gyr {thick lines) and t = 15 
Gyr {thin lines): Ctherm {solid lines), Cmag {dotted lines), and eturb {dash-dotted lines), (b) Variation of 
magnetic energy density with gas density at t = 12 Gyr {dashed line) and < = 15 Gyr {solid line). The flow 
parameters are: Ct = 0.5, SNu = 0.066, Gee = 3 x 10^° and zum = 0. 
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